Coulomb charging energy for arbitrary tunneling strength 
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The Coulomb energy of a small metallic island coupled to an electrode by a tunnel junction is investigated. We employ Monte 
Carlo simulations to determine the effective charging energy for arbitrary tunneling strength. For small tunneling conductance, 
the data agree with analytical results based on a perturbative treatment of electron tunneling, while for very strong tunneling 
recent semiclassical results for large conductance are approached. The data allow for an identification of the range of validity 
of various analytical predictions. 



PACS numbers: 73.40. Gk, 73.40. Rw, 74.50. +r 

Charging effects in metallic nanostructures 0] continue 
to be studied extensively both theoretically and experi- 
mentally. While for the case of weak tunneling the un- 
derlying physics is well understood |^ and starts to find 
mctrological applications j^], for the case of strong tun- 
neling there are several conflicting analytical predictions 
l^-p^. To settle these discrepancies, we have performed 
precise Monte Carlo calculations for a large range of the 
tunneling strength. 

The simplest device displaying Coulomb charging ef- 
fects is the single electron box (SEB) which is formed by 
a metallic island between a tunnel junction and a gate 
capacitor, see Fig. |l|. The external electrodes are biased 
by a voltage source. If the tunneling resistance Rt is 
large compared to Rk = h/e'^ ^ 25.8 kf2, the charge q 
on the metallic island is quantized, q = —ne, where n is 
the number of excess electrons on the island. A simple 
electrostatic calculation for vanishing electron tunneling 
gives for the ground state an electron number n which 
is the integer closest to ricx = CgU/e. Hence, as a func- 
tion of the applied voltage, n displays the well-known 
Coulomb staircase. 

At finite temperature the Coulomb staircase is rounded 
by thermal fluctuations of the island charge. This ther- 
mal smearing is easily understood and seen experimen- 
tally However, there is also a quantum-mechanical 
rounding of the steps, since the island charge is not 
truly quantized. Electron tunneling leads to hybridiza- 
tion of the states in the lead and island electrodes. For 
very strong tunneling, we expect a complete washout of 
charging effects, and the average electron number (n) 
in the box becomes proportional to the applied voltage, 
(n) = riox- As indicated in Fig. |l|, there is a grad- 
ual breakdown of the Coulomb blockade behavior as the 
dimensionless tunneling conductance at = Rx/Rt in- 
creases. 

The classical single electron charging energy of the 
SEB is given by Ec = e^/2C, where the island capac- 
itance C is the sum of the capacitance Ct of the tun- 
nel junction and the gate capacitance Cg (see Fig. ||). 
This quantity can be extracted from measurements in 



the classical regime The strength of the Coulomb 

blockade effect may be described in terms of an effective 
charging energy E* which coincides with E,, in the limit 
of small quantum fluctuations of the charge, at <C 1, and 
vanishes in the limit of strong tunneling, at ^ 1 . Since 
in the latter case d{n) /ducx ^ 1, we may define E* from 
the slope of the staircase at ricx = by 

e: - i?, (1 - a(7i)/anex|„..=o) . 

A completely equivalent definition of E* in terms of the 
free energy F(nox, at) of the SEB reads 
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Several previous articles have made analytical predic- 
tions on E*/Ec as a function of at- For small at, pertur- 
bation theory in the tunneling term may be employed. 
One finds at zero temperature [p^ 

E*jEc = 1 - (l/7r2)at 4- daa? -f ©(a?) . 

The linear term is readily evaluated . For higher-order 
corrections one must take into account that the charging 
energy in the unperturbed Hamiltonian leads to a cor- 
relation of the Fermi liquids in the lead and island elec- 
trodes. Based on the non-crossing approximation (NCA) , 
Golubev and Zaikin § find (see Eq. (8) in Ref. §) 
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where Li2(a;) — — dz ln(l — z)/z is the dilogarithm 
function. On the other hand, from the systematic evalu- 
ation of all diagrams, Grabert MM obtains (see Eq. (2) 
in Ref. § or Eq. (100) in Ref. |ig) 
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The difference between Eqs. (||) and (^) indicates that 
the NCA becomes a poor approximation for large at ■ 
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FIG. 1. a) The circuit diagram of the single-electron box, 
consisting of a tunnel junction in series with a capacitor Cg . 
The junction resistance is Kt and the junction capacitance is 
Ct. b) The number n of electrons in the box as a function 
of the applied voltage U = {e/Cg)nc^ in the Coulomb block- 
ade regime, Rt ^ Rk (solid line) and for strong tunneling, 
Rt <C Rk (dashed line). 

The case of strong tunneling, a* ^ 1, was first stud- 
ied by Panyukov and Zaikin ||^. Based on an instanton 
approach these authors find (see Eqs. (7) and (8) in 
Ref. §) 



E*jEc = a^t exp(-at/2 -t- 7) 



(4) 



where 7 — 0.577 ... is Euler's constant. Very recently, 
Wang and Grabcrt have evaluated the low-temperature 
partition function of the SEB by related path-integral 
methods [Q . For large at , the path integral is dominated 
by niulti-sluggon trajectories whose contribution can be 
summed analytically yielding (Eq. (10) in Ref. [Q) 



E:/E, = 2a3e-"*/2[i + O(hi(a0M)] 
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The exponential dependence of E* on at in Eqs. (^) and 
(^ is in accordance with a recent renormalization-group 
analysis However, the pre-exponcntial factors differ 
by orders of magnitude in the strong-tunneling regime. 
Available renormalization-group techniques do not re- 
solve this problem. 

To clarify the discrepancies and identify reliable meth- 
ods for further analytical work, we have performed Monte 
Carlo (MC) simulations. We start from the path integral 



for the partition function of the SEB which may be writ- 
ten as 



(6) 



where P = l/fc^T is the inverse temperature, and the 
integral is over all paths of the phase (p{t) in the interval 
-/3/2 < T < /3/2 with (p(/3/2) = <^(-/3/2) modulus 27r. 
The action 

contains two parts describing charging of the island and 
tunneling across the junction. The effect of the Coulomb 
energy is contained in 
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(ip + 2inc^Ecf + EcU, 



Sc[ip] - / dr 

J-p/2 L4^c 

The second part of the action jlj] 

/•/3/2 fP/2 

St[(p]^2 dr dT'a(T-T') sin^ 

J-p/2 J-p/2 

describes tunneling. The Fourier transform of ^(t) reads 
ai = —at\uJi\/^TT for ^ D, where uji = 27rZ//3 are 
the Matsubara frequencies and D is the electronic band- 
width. 

By using the definition (|l|) for the effective charging 
energy, we can relate E* to the partition function, 
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To proceed, we write Eq. (|^) in the equivalent form 

00 p 
Z{n,^,at)= e^™- / Z?[,9]exp(-^„[i?]) , (8) 

m— — C30 

where we have separated the functional integration into 
a sum over winding numbers m and an integration over 
all paths 'i?(r) which are now subject to the boundary 
conditions 'd{-(3/2) = 'd{(3/2) = 0. Here, the action 
Sm [t9] is given by 
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^From Eqs. (0) and (|8|), a convenient starting point for 
MC simulation is given by 



K _ E„^^'/^[^]exp(-5„ 



Ec f3E, E™/^['9]exp(-5,„[i?]) 
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which expresses E* in terms of the mean squared winding 
number. 
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FIG. 2. Typical paths fir) generated in the MC simu- 
lation for large at and low temperatures. Specifically, this 
configuration was present after 50,000 MC steps for at — 30 
and f3Ec = 500, taking ip{T) = initially. 

For the MC simulation, we discretize imaginary time 
in P slices of sufficiently small length At = P/P, and 
then sample i9(Ti) and the winding number m from the 
action Sm [&] specified in Eq. (^) . Measuring along the 



MC trajectory allows to extract E* according to Eq. ( 10 ) 
As is apparent from Eq. (^), the MC weight function is 
not necessarily positive definite for non-integer riex- The 
resulting interference can lead to numerical instabilities, 
especially near half-integer values of Ucx- The MC algo- 
rithm developed in Refs. deals with the case of 
general ricx and hence suffers from this problem. How- 
ever, since we focus on the effective charging energy only 
and hence on the value riex = 0, there is no instability 
problem in our algorithm. This circumstance allowed us 
to reach temperatures of about one order of magnitude 
lower than studied in Ref. typically (3Ec = 500. 

For a given parameter set (at,(3Ec), we have first de- 
termined the Trotter number P by empirically checking 
convergence to the large-P limit. Typically, a value of 
P — 5f3Ec was sufficient. Special care is then neces- 
sary for large at since the acceptance rates can be very 
low. The most important paths fir) encountered in this 
regime are of the form shown in Fig. |2[ These paths 
closely resemble the multi-sluggon trajectories contribut- 
ing to the semiclassical result . The predominant 
occurrence of these paths in the MC sampling gives a 
strong indication that the semiclassical sluggon calculus 
is indeed appropriate for large at- The MC updating 
employs single-particle moves, where different samples 
are separated by five passes. Results reported below are 
obtained from several million samples per parameter set. 
The simulations were carried out on IBM RISC 6000/590 



workstations. 

Our data for the effective charging energy are shown 
in Fig . p| together with the various analytical predic- 
tions [|6|,^[- fiOHl4| ] and the MC results by Falci et al. 
p2[ . For weak tunneling, the zero-temperature limit 
was approached already for rather high temperatures, 
f3Ec = 100. However, it turns out to be quite a demand- 
ing task to reach this limit for large at- In fact, the data 
points for at — 20 and at = 25 (which were obtained for 
/3Ec = 500) are not in the true zero-temperature limit 
yet. This follows from comparing additional MC data 
obtained at different temperatures (not shown here) and 
has to be taken into account when comparing the finite- 
temperature MC data with the zero-temperature analyt- 
ical predictions (^) and (||). 



H 




20 25 30 



FIG. 3. MC data for as a function of at and com- 
parison to previous work. The perturbative results for small 
at are depicted as the dashed-dotted curve (see Golubev and 
Zaikin |^) and the solid curve (see Grabert respec- 
tively. The analytical results for large at are the dotted curve 
(for zero temperature, see Panyukov and Zaikin |^) and the 
dashed curve (for (3Ec — 500, see Wang and Grabert 
respectively. The MC data by Falci et al. are shown as 
open circles, and our MC data are given by filled diamonds. 
In the inset, the data for large at together with the curves 
from Refs. [g,^ are replotted on a semi-logarithmic scale. 
The data points were obtained at /3Ec = 100 for at < 15 and 
at PEc = 500 for at > 15. Statistical errors on our MC data 
are smaller than the symbol size unless indicated by vertical 
bars. 

For the case of weak tunneling, the theoretical predic- 
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tions by Grabert [^,0 given in Eq. (^) are confirmed 
with very good precision up to surpringly large values of 
the tunneling conductance, at < 10, while the NCA re- 
sult (H) becomes inaccurate for at > 5. We suspect that 
the discrepancy between the MC data of Ref. Q and 
our results is due to the lower temperatures employed 
here. For at < 15, our simulations have converged to the 
zero-temperature limit. 

For the case of strong tunneling, the predictions of Pa- 
nyukov and Zaikin Q and of the sluggon calculus dif- 
fer by about one order of magnitude due to different pre- 
expontial factors, while our numerical data for f3Ec = 500 
lie in between both predictions ||l9[| . Since a further de- 
crease of temperature will suppress thermal fluctuations 
even more, our MC results give a lower bound to the 
true zero-temperature result for E*. However, the values 
of at are still not large enough to confirm the sluggon 
prediction (^) unambiguously. 

In conclusion, we have developed and carried out 
Monte Carlo simulations for the effective charging energy 
of the single electron box. The data describe a grad- 
ual smooth crossover between the analytically accessible 
limits of weak tunneling (where perturbation theory in 
the tunneling conductance applies) and strong tunneling 
(where a semiclassical calculus is appropriate). 
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